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THEORY OF TWISTED NONUNIFORMLY HEATED BARS 

B.F. Shorr 

An approximate theory of twisted, nonuniformly heated bars of /l^l* 

arbitary cross section is discussed, with the nonlinear distribution 

of normal stresses taken into account, as applicable, according to 

2 

experimental data, up to twist parameter g %5. 

The theory of twisted bars* which is of great importance in calcu- 
lation of airscrews and compressor and turbine blades, has been developed 
in two directions. In the studies of P.M. Riz [1, 2], A. I. Lur'ye and 
G.Yu. Dzhanelidze [3] and others, the problem was solved by the methods 
of elasticity theory. Because of the complexity of the solution, they 
were reduced to final form, only for bars of the simplest cross sections. 

V.P. Vetchinkin and N.N. Polyakhov [4] and I. A. Birger and the author 
(in 195*0 proposed approximate methods, which are applicable to uniform- 
ly heated bar® with specific types of cross sections. A general theory, 
which imposes no restrictions on cross section shape and is valid for 
both slightly and moderately twisted bar s (a refined classification of 
bar s by degree of twist will be established below), with nonuniform heat- 
and with account taken of variable elasticity parameters, is presented 
below. 


1. We consider a bar of constant cross section with, in the un- 
stressed state, a uniform twist relative to the rectilinear z^ axis, 
which passes through a certain point i normal to the cross section axis 
(Fig. 1). 

We limit ourselves to analysis of bars, for which, before and after 
deformation, the following condition holds true 

(«/?)* ( 1 . 1 ) 


“Numbers in the margin indicate pagination in the foreign text. 
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where « is the relative angle of twist of the ^art, 
R is the distance from the z^ axis to the most 
remote point of the section. 


Since 


#r «s» tgt r 


( 1 . 2 ) 


. \ ^ where $ r is the angle between the axis and the 

* screw line which connects the corresponding points 
of adjacent sections, condition (1.1) is equivalent 
* ? * s * 1 to the assumption that this angle is small, which 

permits it to be considered that 


sin8 r »«,ssar, cos8,«i (0 <r</q 


( 1 . 3 ) 


Together with plane of the cross section n n » we introduce "orthogonal 

1 0 
section" n> a surface, the shape of which would be the cross section, 

if the bar were reduced to the twisted state (o»o Q ) **rom the untwisted 

state (ot*0) by free twisting (Pig. 2) . 



By definition, the orthogonal section 
coincides with the warp surface at a*aQ and, 
with condition (1.1), it is described by the 
equation [5, 6] 


in which 


u\> =*«..?(*. y) (l* 1 *) 

(«<.?*')* < 1. (*o?»')* < 1 (1.5) 


Pig. 2, 


I f 

where <|>(x, y) is the twisting function, $ *3<j>/3x,4 »3<f>/3y; x, y, z 

A V 

is the rectangular coordinate system which moves together with the cross 
section. For bar s with variable elasticity parameters E(x, y) and 


The term "orthogonal section" was proposed by I. A. Birger, as applied 
to bars of elongated cross section, for which section II can be identified 
approximately with the surface normal to the screw lines. 
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( 1 . 6 ) 


and G(x,y), function <f> should satisfy the equation [7] 

(o?,v + “ G * (y — y*) - <V (* - **) 

and the boundary conditions on the profile 

(?*' — y + yi )dy — (?*' + * — *0 dx =* o (1.7) 

Because of condition (1.5), no distinction can be made between the 
areas and equations of the profiles of the orthogonal and plane sections. 
To be definite, we assume that a given orthogonal section n corresponds 
to the plane section n o , with respect to which 

^OfdP-O a. 8) 

Eq. (1.4) characterizes the configuration of a twisted bar, to 
within the position of twist axis z^. The final results do not depend 
on the direction which satisfies conditions (1.1), selected as the 
direction of the axisj however, as for prismatic bars, the theory is 
simplified if it is assumed that the z^ axis passes through the center 
of rigidity of the section. 

Since, corresponding to twist a , the initial "shifts" 

T«° - * 0 (<P*' - y + yi), T«w° = *0 (?/ + * — *<) (1.9) 

should not depend on selection of the position of the axis, 

— Ifi, <Pv' «■*„' + *, (1.10) 

Here, $ is the twisting function, determined as usual, on condi- 
tion that the axis of rotation passes through the coordinate origin. 

2. With condition (1.1), the elastic twist deformation can be con- 
sidered as a continuation of the initial "deformation," and the shifts 
can be determined by the formulas 
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fj *= — 0 (*/ — !/t). u 1 = O9 (2.1) 

where 0 is the relative angle of elastic twisting which, at small 
Initial twist angles, can be commensurable with a Q . 

The twisting stress In an orthogonal section is 

x ** w y + Vi), =* co (<f v ' .f. x — *i) (2.2) 

For determination of the normal stresses, we consider that, at 
all points, the profile of the orthogonal section is normal to the 
open lateral surface of the bar. Therefore, it can be proposed that 
the tangential stresses on this section, which are not connected with 
twisting (at 0*0), are of secondary importance (just as for a plane 
section in the theory of transverse bending of beams). Therefore, as 
the initial direction of a longitudinal fiber of a twisted rod N Q 
(Fig. 2), it is natural to assume the direction normal to the orthogonal 
surface with the directing cosines 

U “ cos (N 0 *) » — a„ 9 x', m 9 - cos (N Q y) = — ot a <f v ’, n 6 - cos {N Q z) =1 (2.3) 

and, following the general theory of bars, to disregard the pressure of 
one longitudinal fiber on another, i.e., to consider 

<>•*£(• — T «) (2.4) 

where e is the longitudinal stress and yt is the thermal expansion of 
the fiber. 

It is evident that, in the general case, the directions of the 
longitudinal and screw fibers do not coincide. 

The length of a unit fiber in the initial state is 

*0 - * FT+ V + V * dt [1 + i. # 0 * ( 9x <* + 9( /*)] (2.5) 
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( 2 . 6 ) 


and, after deformation, 

+ •.)• + <■+"»* • 

where e z is the deformation of the fiber parallel to the rod axis, *« 
* 0 +AA, m«m o +Am and A£, Am are the increments of the angles of inclina- 
tion of the longitudinal fiber, in connection with elastic twisting of 
the bar. 

Based on the ’’orthogonal sections’’ hypothesis similar to the 
hypotheses of plane or Irregular sections, we assume that, in stretching 
and bending, the orthogonal section does not change its shape, moving 
according to the laws of solids. Then [8], 

•«"■# — **x — i v y+ 69 (2.7) 

where e Q , < x , are components of the plane of deformation, 8«d0/dz, 
and 

A/- — 9(y — y t ), Am ■ 6 (* — *,) 

( 2 . 8 ) 

In conformance with (2.3) and (2.8), Eq. (2.6) takes the form 


dtssdi{ 1 +«, H- 1 (*, 9 / *0(y — y*)l* + 1 I*#9k' — ®(ar — x,)]*} (2.9) 

from which 

LOV (,2.10) 


where polar coordinates r, ip, with the pole in the center of rigidity, 
are introduced, and it is taken into consideration that 

to- ?*' (2.11) 

The normal stress in the orthogonal section with direction N(i,m,n) 
is 


• »£(«,_ M 4. fcp - 1 0 V* - ft) 


( 2 . 12 ) 
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Besides torsional stress tangential shear stresses T ^ act 
in the orthogonal sections. The latter are connected with stress o by 
the equation of equilibrium 


*• , , V <W 11*/! + ,, „ o 

* .in * • 


- 4 - 

d*T 9 * 


Oy 


(2.13) 


where p is the bulk force component in the direction of the longitudinal 
fiber. 

The total tangential stresses 

(1) (2) 

T "T +T 

The stresses in the plane of the section can be found from the 
conditions of equilibrium of a longitudinal element bounded by sections 
II and n Q (Fig. 3), which gives 

Ipx — r)j — ,"4 8 »] 0 - *' 1 « „"'i ‘{(i/J — fi) 0 + o — 

(0 + •»-») ( 2 . 14 ) 

+ /(*•»*)! » + +/4« 4 ) » + o » ,0 



Fig. 3- 


In particular, in the case of stretching 
of the bar, by dropping out all nonlinear terms, 
we obtain, with x^*y^“0 and t»0, 

*" “**•?• + if 

«.«[-*+£ (»,'+.)] (2 ‘ 15) 




Here, I p is the polar moment of inertia, and 
T is the geometric rigidity in twisting, which 
coincides with the precise solution of P.M. Riz 
[ 1 ]. 


At points where <j>"0, Eq. (2.14) coincide with the conventional 
formulas for the stress in inclined areas. 


With the force factors in the plane of the section considered to 
be known, we write the equilibrium conditions in the form 


nM 
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Pam j a ' dF ' ^ t.,W, Qy - J T ,SdF 

*/-- A/v — | xMF, M ti - ^ [r,„«(x -*,) - z\ x (,j -' yi )\dF (2.1 6) 

We note that, between the force factors In a twisted bar* there are 
the relationships 

$«- ~~ar+ 

0,-^* -«(/>**+ A/,) (2.17) 


which Eq. (2.16) identically satisfy. 

By substituting the values of the stresses in the first and last 
three equations of (2.16), we obtain a system of four equations, which 
are linear with respect to the components of the plane of deformation 
e Q , < x , Ky, and nonlinear relative to angle 6. Without presenting the 
general expressions, because of their bulk, we consider some applica- 
tions of theoretical or practical value. 

3. For profiles of dimensions x, y of one order, tangential 
( 2 ) 

stresses x' ' and the constrained twisting effect can be disregarded. 
In this case, we arrive at the following system of equations: 

. /> + />,- E m !*,/•’ - * X S V - H-0[« 9 (/„ - T) + [ 0 /„ ]} 

"f" ^*1 xv'~*yl 0 (/ px — T t ) O/pijl 

My *f* Myt ■ — E„ | — X«/y Y.yJ 8 (I py T y) -{- ' 0/ 

M, - G m n + E m {«, [« # (I p 0/ P l - X, [« # (Ipy - Ty) + 0/pyJ- 

— *¥ (®» (Ip* — T*) + ®7p*J + 0a o [# # (/, — T r)-f- jj(/ r — 7%) J + J — 

— «,(/?< — D t y) — QBt 


(3.1) 


Here, we Introduce designations for the reduced geometric 
characteristics 


Si - ^ydF, - ^EPxdP, 

■ /« - /,- ^ITz'dP, tn-^EPxydP 
Ii-^E'r'dP, /„. JjEVtyf/?, /„- ^E*r'xdP, I,-^E*r*dP 
T - rW IT. - 'W T v ~ ^E*xM+r*)dP 

TV- T„.^Fr»i f t;+ r ')dF 
^ j ^EdP * 2 (1 + (*) G,*, court j 


(3.2) 


and temperature factors 


P, - ^E-fldF, a /„- Y'^ dF ' " “ ^E-itxdP 
Bt -■ Z? (f ■ fa’ + r s ) rfF 


(3.3) 


In the derivation of system (3-D* it was taken into account that, 

because of the properties of twisting function $ and the assumption 

w Q «R, the role of the terms which reflect the effect of tangential 

stresses on quantities P, M , M is insignificant, within the limits of 

x y 

accuracy of the theory. For simplicity in recording, subscript i in 
moment M x is dropped subsequently. 


By proper selection of the coordinate origin and directions of the 

x, y axes, condition S “S “I “0 can be ensured. By introducing the 

a y xy 

principal values of the components of the plane of deformation 

l 

••-«.+ *[«. ('p-n+Jo/p] 

*«•-*«- ^[«,(/pp-r,)+ Jo/, th ] (3.*0 

*V) + !•/,„] 


<**!£*** 
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the principal polar coordinate 


rS-r'-'f-'-fXx-'jfy (r%5 0) 


(3.5) 


'1H6 


and the principal value of function 


/ T I _ T I T 

(9 V - + *— + -X-jf- x •(- -*Tpl y 


(3.6) 


which satisfy the conditions 


I »mi *-1 *-0. 

P px py * 


T»«T *"T *"0 

x y 


where I *, T * (n Is an arbitrary subscript), are determined by F,q. (3*2), 
n n o o 

with substitution of r# for v c and ($*)^’ for we represent the 

first three equations of (3*1) in the form 


- p + p, „ 


-BJT 


(3.7) 


Eq. (3*7) are the same in form as the corresponding relationships 

for rectangular bars. However, in twisted bars* parameters < * and ie * 

a y 

do not coincide with the components of the elastic curvature of the 
rod axis or e* with its elongation. 

It follows from (2.12), (3.*0-(3.7) that, with 0»O, 


f P + Pi W-i 4 1 

£ Ix 7 — t'J 


(3.8) 


The last equation of (3*1) can be presented in the form 


•» — • 4 (1 + v + P*) + 3 (1 + |i) *,0* — — y — -r (1 + 1 4 ) ®* “y - 


(3.9) 
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where 


ft - *• 
8 »" ITT 


v* _ 


7 H “ 7^ ( ; w - 7*) + ^ (/*, - 7V) - (/?,* - 




CT 


/* r K » 

(r 7 *- 77^+7: 




( 3 . 10 ) 

(3.11) 

(3.12) 


where B * and B * are determined by Eq. (3.3), with the substitution 
2 t 2 t0 

of r* for r and (<f>*)^' for and 


P*- 2(1 -MO 


(3.13) 


With a Q - 0 (Initially untwisted bar), Eq. (3*9) changes to 


/.• 


(i + f*)-y-0* + (| + v)0-,0. 


(3-1^) 


For cases of stretching and twisting bars of circular and elongated 
rectangular cress sections, with t-0, formulas are obtained from (3*1*0, 
which coincide with the results of S.P. Timoshenko [9]. With© o *0, 

Eq. (3*1*0, besides the trivial solution 0*0, has a second solution 


•-± 


/ 


7/<i +c) 


(3.15) 


which takes on a real value at \j<-1. /1^7 

By assuming in (3.12), the value v"-l, we find the condition 
for loss of stability of a rectangular bar due to twisting. At M x « 

My»0, t«0, this gives the known formula for the value of the critical 
force [10] 
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p « GTF 

P * V 


2 3 

By disregarding the 8 and 0 terms in Eq. (3»9)» we solve it 


for 0, 

• ~ k t ( 0 , - ( *. - 


(^-1) 


The formula for the effective geometric torsional rigidity of the 
twisted bar results from relationship (4.1) 

T 0 «T(l+v+$ 2 ) 

which, at v<<l+g? , coincides with the formula obtained for this case by 
Chen Chu [11], as applied to hollow tubes with an elongated blsymmetrlcal 
section. 

2 

Dimensionless coefficient g plays an important part in the theory 
of twisted bars, and it car. be called the twist parameter. Depending on 

p 

the value of g , it is expedient to provisionally divide twisted bars 
into three groups. 

2 

a. Slightly twisted bars, g «1. In this case, with 6 Q -0, elastic 

twisting angle 6 is proportional to angle a., strains e and e_ practi- 

o z 

cally coincide, and normal stresses a can be calculated by the conven- 
tional formula for a rectangular bar. 


For a section with two axes of symmetry, with t»0 and v<<i(k 0 ^l), 
we have 


%m TBT *)• ^ “ P77 


(4.2) 


which coincides with the results of G.Yu. Dzhanelidze [12] . 


2 

b. Moderately twisted bar 8 * she value of 0 is commensurable 
with unity. With Increase of g 2 , the torsional rigidity of the bar 
increases considerably, which is confirmed by experiments [11], the 
tensile and bending rigidities decrease (for asymmetrical profiles). 
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and the normal stresses In the orthogonal section (and the stresses 
practically equal to them In the plane of the cross section) are 
redistributed according to relationship ( 3 * 8 ). 


As was shown by the tensile tests of twisted bars by the author 
[ 13 ] * the approximate theory reported, based on assumptions ( 1 . 1 ), 
gives good accuracy up to 02 ^ 5 , 

2 

For blades and air screws, the value of $ can reach 2-3 or more. 


Together with the concise notation of (3.7) and (4.1), the basic 

system of equations of a twisted bar of arbltary cross section, with 
2 2 

0 45 and v«l+g , can be represented In the following developed form 

4 

*i m 2 (4.3) 

*« i 


where Cj should be understood to be strain components e Q , te x , < y , 6 and 
L^, the force and temperature factors P+P t , M y +M yt , M x +M xt* /148 

V a o (B t- B t^» 40(3 c0 ® fflcients a jj c » which satisfy the theorem of 
reciprocity, have the following values 

(4.4) 




(/ .. }‘)M / t 

I l-?(l i »•)**«** • Vf • 


m y 






tL 


\ 


V P»* *1 


/ — T 

*9M 


•m | • t ^(H !*)*,«.* *** *■ — *** *■ ■* 


2 

c. Strongly twisted bars, 0*5. A general theory of such bars 
has not been developed. 
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It follows from Eq. (3.2) that. In the case of slightly warped 
profiles, for which >^0, the equality T n ;vl n occurs, and all the terms 
which contain angle a Q revert , to zero, i.e., in this case, the initial 
twist does not play a part. An example of such a bar is a round 
cylinder. 


5. Twisting has the greatest effect on bars with elongated, strong- 
ly warped profiles, for which the value of integrals (3.2) depends 
basically on the values of the subintegral functions at points remote 
from the rod axis, where r»h (h is the greatest thickness of the pro- 
file), and the initial "shifts” in the plane normal to the radius 

Ti+° = ®o ( r + t?*) (5.1) 

become small, compared with the initial "bends" of the fiber 

“W = Y*o( r ~ T?*) (5.2) 

With the assumption that Y z ^°^» we (Fig* 2) 

<p*' =s — cos (A' o 0 - — y- ?*' ~ v = $<• (5*3) 

i.e., at points sufficiently remote from the axis, the longitudinal 
fiber coincides with the screw line, and the orthogonal section, with 
a section of the corresponding screw surface, as was adopted by I. A. 
Birger. 

Because of equality (5*3), we have 

?’n ~ 0, (?),; T n • « o, « o, u t ; ^ o, v* « v 

and all formulas of Sections 3 and 4 are significantly simplified. 

Instead of relationship (2.3), there will be 

4,^ «„(.'/ 'J>)< m «*~ a ^ 
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and, Instead of (2.12), with ©<<0^, 


o : A* (5„ v x T ■ *„;/ I - (5s | xjir' — f/ 


(5.5) 


If, to condition (1.8) adopted above 

| /%,//-• ... o \ Eyrdb ' *, 0 


(5.6) 


is added, formulas (3*7) remain in force, and conditions' (5.6) are 
used to determine the coordinates of the center of rigidity. 


JEJy the substitution of (5*5) in torque Eq. (2.16), without con- 


2 3 

sideration of the 0 and 0 terms, we obtain 


, 2 (Hu) p 

0 o = =*o v = °( I *f- v + p 5 ) -|- G~f ■[• f « o 0 ) 7ifr*dP 


(5.7) 


where M ^ is the torque of tangential stresses t^, which can be 
z 

found for elongated sections without determining the stresses them- 
selves, if it is assumed that 6 1 ccr” 3- * ' , i.e., that warping is re- 

r T m 

duced primarily to rotation of individual elements of the cross 
section relative to the radius. Then, 


=» — (* — Xi ) 

Mz'" = - 5 l(t„«»?) f '+(x IV <*»9)w') dF + \ 9((t„ <*>)*' + dF 

F F 

« - § ?[^x {t) dy - t f ,<s>rf*l - \ ? d J. dF - \ ?P dF 


(5.8) 


Because of the boundary conditions on the profile, the first 
integral equals zero, so that, with 9t/3z«0, we have 


M, 


.«’> - - FmUO - 5 E*?r*dF - $ 9 pdF 


(5.9) 


where 
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^ E°fdF 


(5-10) 


By substitution of (.5*9) in (5.7), we arriye at the differential 
equation of constrained twisting 


where 


0 X*0 = — Ar # X* (0 o — a 0 v C) 


(5.11) 


x 5 r... 

+ tov„ ’ 



(5.12) 


with the boundary conditions 0=0 in the end connection and 0=0 at the 
free end. 

With o q = 0 , p=0, t=0, Eq. (5.12) coincides with the equation of 
constrained twisting of the theory of thin bars [8] . 

By designating 

= - EmlX M, « , A// = M t - M, -f 5 - C m 7«„v (5-13) 

F 

we obtain jif • 

• = (5.14) 

m 


and, consequently. 


/150 


/■•! /■, ^ t M x -|- Al xl 


It, 


x -i — 


y — 7 ~ ? + 


2 ( H- n) A 0 * # r* — j 


For thin wall profiles, on the basis of (5-8), 


«= (y — y*) ~(x — x i) y,' 


(5.15) 

(5.16) 


(5.17) 


where s is a coordinate read along the midline of the profile. 
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After integrating (5.17), with changes of $ through the profile 
disregarded, we obtain, to within a constant. 


— w 


(5.18) 


where w is the doubled sector area. 

to) 

Tangential stress x , which acts along the midline of the profile, 
is 


.<*> 


7^*) [\ P ('») h (*») *1 + 


K-j •(»*)* (0*»1 

0 J 


(5.19) 


which, with 3t/3z=0 and with (2.17) taken into account, gives 


where 


t„<*> » _L /£ifi 


mi ^ r 1 S *«*«*-? p(s l )h(M^dt l + 

A _ * «* 


+ 


V*w . 


K 
l J 


, [</■* -M.) _ (Pl , +A/ >) £^> + B _ ^1]} 


F (*) “ S *y, (*) « \ E°/iyds lt Sy ( s ) = ( 

? • ■ i 

i • 

(*) = J FWf,, 7 P *(i) =0 ( E-hr'dSi 

' a y 


(5.20) 


(5.21) 


and S q is the length of the midline. 

6. In a number of cases, the constraint effect appears only near 
fastened cross sections, where 6=0 and the initial twisting has no 
effect but, in sufficiently remote sections from the end connection, 
the twisting can be considered practically free, which considerably 
simplifies the calculations. 

Without consideration of constraint, Eq. (3.9), for elongated 
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( 6 . 1 ) 


profiles, Is reduced to the form 

p + 3?S» + 2(l + v + fi s ){.2($ 0 -pv) So-^0,) 


Analysis of Eq. (6.1) which reduces arbitrary cases of deformation 

of twisted bars to the single relationship £=£(£ v,f3), confirms that 

5 ■? 0 

the role of nonlinear terms £ , £ is small, as a rule. 


By using the equations of elasticity theory i.n a nonorthogonal 
curvilinear coordinate system, V.M. Marchenko developed, in general 
form, a theory of free torsion and tension without bending of a uniform- 
ly heated twisted bar, with an arbitrary value of angle a Q , in which, 
for an elongated elliptical section with semiaxes ab(a>b) andy=l/3, 
the solution in approximate form (by the method of variations) was 
reduced to calculation formulas, which have the form (in our notation) 


0 




( Pa t \ 

Vcf ~ laj-n) 



i + 


5 a 
28 


ifT 


*0 ~*l 





( 6 . 2 ) 


For 

formulas 


an elongated ellipse (n<l), with H=l/3 and v<< 1+g 2 , 
(3*13), (4-1) and (4.4), we have 


k „ _ *• (a, 31«*ao'\ 

* 0 L 1 45 i) J • p'to •») J’ 




*0 

CT 


®u - ®4i = — *• 


by 


and comparing solutions (4.3) and (6.2) shows that they practically 
coincide. The normal stresses in the center of the section are: 


according to V.M. Marchenko 


o (0) -*,[(« + 


60 o»a 0 *\ P 29 M,-\ 

224 yj JF 58 *° Fi\ J 


according to formula (5.15) 


(6.3) 


/151 
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The divergence of the results also is small. 


Thus, comparison of the approximate theory with experimental data 
and with the partial solutions obtained with the use of the system of 
elasticity theory, indicates its sufficient generality and accuracy. 
The author thanks I. A. Birger for discussion of the work. 
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